Evidence for a large magnetic heat current in insulating cuprates 
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The in-plane thermal conductivity k of the two-dimensional antiferromagnetic monolayer cuprate 
Sr2Cu02Cl2 is studied. Analysis of the unusual temperature dependence of k reveals that at low 
temperatures the heat is carried by phonons, whereas at high temperatures magnetic excitations 
contribute significantly. Comparison with other insulating layered cuprates suggests that a large 
magnetic contribution to the thermal conductivity is an intrinsic property of these materials. 



There is growing experimental evidence that spin ex- 
citations may contribute significantly to the heat cur- 
rent in low-dimensional spin systems. This seems to be 
well established for one-dimensional (ID) systemsi*2iMi^. 
For example, in the insulating spin-ladder material 
Sr 1 4_ 2 .Ca 2; Cu2404i a large magnetic contribution k m 
to the thermal conductivity k can be derived from a 
pronounced double-peak structure of k along the lad- 
der direction^. The situation is less clear in two- 
dimensional (2D) spin systems. These are, however, of 
particular importance due to their relevance for high- 
temperature superconductivity^. A double-peak struc- 
ture comparable to that in ID systems is found in the 
in-plane thermal conductivity of insulating 2D cuprates 
such as La 2 Cu0 4 (LCO) and YBa 2 Cu 3 6 (YBCO) (Ref. 
Iil l9ll0j) . This may indicate a sizable magnetic contribu- 
tion to the heat current at high temperatures-:. How- 
ever, the phononic thermal conductivity k p h may show 
a double-peak structure also, as a result of pronounced 
(resonant) scattering in a narrow temperature range. 
Such scattering may arise from the presence of local mag- 
netic excitations, as was recently shown for the 2D spin- 
dimer system SrCu2(B03)2 (Ref. Hl|) . or it may arise 
from the presence of soft phonon modest. The latter 
was suggested for LCO and YBCO, in which soft modes 
e.g. associated with tilt distortions of the CuO polyhedra 
are known to be presen1j2*i£. An additional complication 
arises from a strong sensitivity of the double-peak struc- 
ture to light oxygen doping 9 . 

A material of particular interest in this context is 
Sr2Cu02Ci2 (SCOC). It is structurally very similar to 
LCO: It contains Cu02-layers as in LCO, but the out-of- 
plane oxygen ions at the apices of the CuOe-octahedra 
are replaced by CI and La by Sr. The material has several 
advantages compared to LCO and YBCO (see e.g. Ref. 
6): (1) SCOC does not exhibit any distortion from tetrag- 
onal symmetry down to at least 10K so that there is no 
structural instability associated with soft tilting modes. 
(2) Because of the absence of tilt distortions the magnetic 
properties are simpler than those of LCO. For example, 
there is no Dzyaloshinski-Moriya exchange interaction. 
Thus, SCOC is believed to represent the best realization 
of a two-dimensional square-lattice S = 1/2 Heisenberg 
antiferromagnet. (3) In contrast to LCO and YBCO, 
SCOC cannot be doped easily with charge carriers. 



In this paper we present measurements of the in-plane 
thermal conductivity k of SCOC. We identify a double- 
peak structure from a pronounced high-temperature 
shoulder around 230 K. Analysis of these data and a 
comparison to LCO and YBCO shows that it is very un- 
likely that the double peak structure arises from anoma- 
lous phonon damping due to scattering on soft lattice 
modes or magnetic excitations. The data indicate instead 
a large magnetic thermal conductivity at high tempera- 
tures as an intrinsic feature of the insulating 2D cuprates. 

We studied a single crystal of Sr 2 Cu02Cl2 of rect- 
angular form (1x3x4 mm 3 ) with the short direc- 
tion along the crystallographic c axis. It was grown by 
the traveling-solvent floating zone method. The thermal 
conductivity was measured with the heat current within 
the Cu02 planes by a conventional steady-state method 
using a differential Chromel-Au+0.07%Fe-thermocouple. 
Typical temperature gradients were of the order of 0.2 K. 
The absolute accuracy of our data is restricted by un- 
certainties in the sample geometry whereas the relative 
accuracy is of the order of a few % 

We show in Fig.^the in-plane thermal conductivity of 
SCOC as a function of temperature. We identify a max- 
imum at « 30 K and a shoulder at high temperatures 
around 230K. The pronounced low-temperature maxi- 
mum of k indicates a high crystal quality. We note that 
k is independent of a magnetic field (< 8 Tesla) applied 
within the Cu02 planes perpendicular to the heat cur- 
rent. For comparison we show in Fig^the in-plane ther- 
mal conductivity of a single crystal of LCO measured by 
Nakamura et al£. These data also reveal a double-peak 
structure. The absolute value of k at the low-temperature 
maximum is smaller than in SCOC. One reason may be 
that LCO is more sensitive to defects, resulting e.g. from 
excess oxygen, which introduces lattice defects and hole 
doping and thus reduces the mean free path of the heat- 
carrying excitations. 

Both compounds, SCOC and LCO, are antiferromag- 
netic insulators. In an insulator the heat is usually car- 
ried by phonons. The typical behavior of k v h of a crys- 
talline insulator is shown by the solid lines in Fig. ^ 
These curves represent fits to the low-temperature max- 
imum of k (fitted below about 50 K) of SCOC and LCO 
using the standard Debye model for the thermal conduc- 
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FIG. 1: In-plane thermal conductivity k(T) of S^CuChCb 
(circles) and La2Cu04 (dotted line; data from Nakamura et 
alp-). Solid lines: fits to k p h using the Debye model—*—. The 
fitting parameters for S^CuO^Cb (I^CuCU) are: _D/lCP 17 s 
= 3.4 (2.6); P/f0~ 43 s 3 = 0.15 (21); [//1(T 3 V/K = i.o (2.2); 
u = 4.9 (4.4). The point defect scattering (P) is smaller in 
SCOC. Inset: k m = k — k v u for S^CuOaCla (circles) and 
La2Cu04 (dotted line) (see text). 

tivity of acoustic phononsii*i^: 

o 

Here Gd is the Debye temperature and v p h the sound 
velocity. Due to the lack of experimental data for SCOC 
we use for both compounds the values reported for LCO 
(Q D « 385 Ki£; v ph w 5.2 • 10 3 m/s 16 ). uj is the phonon 
frequency, x — fku/kBT, and t (x, T) is the phonon re- 
laxation time given by 

r -i = M + D ^ + P ^ + UTcj 3 exp ( . 
L \ uT J 

The four terms refer to the scattering rates for bound- 
ary scattering, scattering on planar defects, on point de- 
fects, and phonon-phonon Umklapp scattering, respec- 
tively. L w 1mm is the sample length, and D, P, U, 
and u are fitting constants. The low-temperature data 
are described very well by these fits. The fit parameters 
are given in the caption of Fig. ^ The decrease of k p h 
at high temperatures is due to phonon-phonon Umklapp 
scattering. 

For several reasons it is very unlikely that the high- 
temperature increase of k is due to conventional heat 
transport by phonons: (1) The contribution to k from 
acoustic phonons, as described above, decreases at high 
temperatures. (2) The contribution of optical phonons 
to the heat current is usually much smaller than that of 
acoustic phonons, even in compounds with a very large 
number of atoms in the unit celi£*i£, so that heat trans- 
port by optical phonons is very unlikely to cause the high- 
temperature maximum. (3) The out-of-plane thermal 



conductivity k c of LCO behaves as the in-plane thermal 
conductivity k at low temperatures, but k c shows no indi- 
cation of a high-temperature maximum^. Such strongly 
temperature dependent anisotropy is not expected for 
purely phononic heat conduction. 

Additional phonon scattering, active in a narrow tem- 
perature range close to the minimum of k, may in 
principle cause a double-peak structure. However, res- 
onant scattering on local magnetic excitations as in 
SrCu2(B03)2 (Ref. lTlj) cannot be the correct explanation 
in the present case: In the 2D square-lattice cuprates the 
dispersion of magnetic excitations ranges from w to 
2J /ks > 2000 K (J is the in-plane exchange constant) 
so that there is no reason that scattering on magnetic 
excitations should be most pronounced in a narrow tem- 
perature interval around 100 K>ii*iii Note, in particular, 
that in the 2D cuprates scattering on magnetic excita- 
tions should not disappear above the Neel-temperature 
Tn, because the relevant energy scale is set by J ^> Tjv 
(see below). Additional phonon damping from scattering 
on soft lattice modes as suggested in Ref.^is also unlikely 
as a cause of the double peak: (1) There are no lattice 
instabilities in SCOC, rendering this mechanism unim- 
portant for this material. (2) The double-peak structure 
is also present in the tetragonal low-temperature phase 
of Eu-doped LCO, in which no soft tilting modes should 
be present either—. (3) The absence of a double-peak 
structure of k c in LCO (Ref. 8) implies that anomalous 
phonon scattering would have to be active only for k. 
Such strong anisotropy of the phonon-phonon scattering 
is not expected. Finally, note that the finding k c < k 
in LCO (Ref.|g) provides evidence against any scattering 
scenario as a cause of the double peak structure: for such 
scattering, if active only for k, but absent for k c , implies 
k < k c , in contradiction to the experimental results. 

The data of Fig. Q] (in particular k > k c ) are most nat- 
urally explained, if an additional channel of heat trans- 
port for the in-plane thermal conductivity is present. 
In an undoped insulating 2D Heisenberg antiferromag- 
net with an electronic gap > 1.5 eV, the only candidate 
for heat transport next to phonons are magnetic exci- 
tations. Their thermal conductivity k m adds to that of 
the phonons, i.e. k = kph + k m . In order to extract k m 
from the data we subtract k p h as obtained from the fit 
of the low-temperature maximum. Note that k m cannot 
be obtained at T < 100K in this way, because Eq.^ was 
fitted to the total k below 50 K. Remarkably, k m is of 
comparable magnitude (roughly of the order 10 W/Km) 
in both compounds (see inset Fig.^). The maximum of 
k m is at w 245 K in SCOC and at » 285 K in LCO. 

Is a magnetic contribution of this size reasonable? We 
estimate k m using the kinetic equation in 2Dia: 

km ~^C-mVm^m ■ (2) 

Here c m is the magnetic specific heat and £ m the mean 
free path of the magnetic excitations. The velocity v m of 
long- wavelength spin waves is wf n coc ~ 1.06 • 10 5 m/s 
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285 
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4.9 



TABLE I: Position Th of the high-temperature maximum of 
k mi Neel-temperature Tjv, in-plane magnetic exchange cou- 
pling constant J, and the ratio J/kaTu for three insulat- 
ing 2D cuprates. Note that J/ksTn is ver y similar in all 
three coin])! nun Is. Tn is from Refs. l(l l7l32L The values of 
J are derived from two-magnon Raman scattering and in- 
frared bimagnon-plus-phonon absorption data*-*—, where the 
Oguchi correction has been taken into account*— The data 
for k in YBCO and LCO are from Refs. |]and|| respectively. 



FIG. 2: Magnetic mean free path l m (T)/a (circles: 
Sr2Cu02Cl2; dotted line: La2Cu04) as extracted from the 
data shown in Fig. (see text), and in-plane magnetic cor- 
relation length £ m (T)/a (solid line: Sr 2 Cu02Ci2; dashed 
line: La2Cu04) as obtained from neutron scattering-. Here, 
a ~ 3.9 Ais the lattice constant. Inset: Specific heat c m /kB 
of a S — 1/2 square- lattice Heisenberg antiferromagnet (see 
appendix) . 



and v m co ~ 1.16 ■ 10 5 m/s, which is obtained from 
v m = V8SZ c Ja/h. Here Z c is the Oguchi correction^ 
Values of J for various 2D cuprates are given in Table 
I. Note that v m is much larger than v p h, as a result of 
J » fegQi). For the specific heat c m we use the theo- 
retical result shown in the inset of Fig. [3 which comes 
from the extrapolation of the high-temperature series 
for the partition sum (see appendix). The maximum is 
c max = 0.4612(5)iVfc B at fc s T max = 0.5956(1)J. Given 
v m , Cm j and using k m as shown in Fig. ^ we obtain an 
estimate of t m (T) using Eq. J3J. £ m decreases strongly 
with increasing temperature (Fig. At room temper- 
ature, £ m /a m 30 for SCOC and « 75 for LCO. These 
values are not excessively large - much larger values of 
£ m have been found in one-dimensional spin systems^*^ 
- rendering a magnetic contribution to the heat current 
in SCOC and LCO very plausible. 

For a better understanding of k m it is instructive to dis- 
cuss the magnetic correlations in the quasi-2D cuprates, 
in particular, the in-plane magnetic correlation length 
£ m (T). In a 2D antiferromagnet long-range order with 
£ m = oo is restricted to T=0. With increasing temper- 
ature spin-flips (or magnons) are excited, which reduce 
£m by breaking the long-range correlation. In the quasi- 
2D materials considered here, the finite magnetic order- 
ing temperature Tn is determined by the mter-plane 
interactionSii, which is much weaker than the in-plane 
exchange interaction J, so that UbTn <C J (see Tab. I). 
In the ordered state at T < Tn, £ to = oo. For T > T N , 
£ m is still large because of the large J. We show £ m (T) 
of SCOC and LCO as inferred from magnetic neutron 
scattering 7 in Fig. Above Tjy, £ m is indeed much 
larger than the lattice constant. However, £ m decreases 
strongly with increasing temperature, approximately ac- 



cording to 7 . £ m (T) ~ exp(27rJ/A: B r). 

From these considerations we may draw several conclu- 
sions for the magnetic heat current: (1) We expect that 
k m is determined by the large in-plane exchange inter- 
action J and not by the inter-plane interactions. There- 
fore, no significant anomaly of k m is expected at Tn- 
(2) At least above Tn the heat-carrying magnetic ex- 
citations are not the familiar collective excitations of a 
magnetically ordered state (i.e. conventional magnons), 
but rather magnetic excitations (of triplet character) in 
a spin-liquid state. Note, however, that also at T < Tn 
the nature of the magnetic excitations of the 2D cuprates 
is under intensive debato 24 i 25 i 26 i 27 ' 28 . (3) The strong, ex- 
ponential decrease of £ m above Tn suggests a similar de- 
crease of £ m (T). This would explain, why the maximum 
of k m occurs at a temperature much lower than that of 
c m : the increase of c m (T < T max ) is overcompensated 
by the strong decrease of £ m (T). Note, however, that the 
results on Sri4_ 2; Ca a ;Cu2404i show that in a spin-liquid 
£ m may be significantly larger than £ m A^ 

A double-peak structure of k has also been reported 
for the insulating bilayer compound YBa2Cu30@- As in 
LCO one observes pronounced anisotropy, i.e. k c does not 
show a high-temperature maximum. Given the existence 
of a large k m in the monolayer cuprates LCO and SCOC, 
the high-temperature maximum of YBCO is also likely 
to be of magnetic origin. This is consistent with the 
systematic variation of the temperature Tjj of the high- 
temperature maximum of k with J for the three different 
insulating cuprates as shown in Tab. [J 

The data in Fig. ^ show that the high-temperature 
maximum is more pronounced in LCO than in SCOC. 
A related observation is the weak high-temperature 
anomaly of k in Pr2Cu04. 29 In view of their rather 
similar magnetic properties one would expect a sim- 
ilar behavior of k m in these compounds. 30 How- 
ever, doping with mobile charge carrier a 9 ' 10 or static 
impuritiesiS^i influences the high-temperature maxi- 
mum strongly. In particular, the high-temperature max- 
imum in YBCO depends in a non-monotonic way on the 
oxygen concentration^. Thus the different behavior of 
k m might result from the fact that LCO and YBCO can 
easily be doped with charge carriers via variation of the 
oxygen content, whereas this is not possible in SCOC. 
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Further experiments with a detailed control of the charge- 
carrier concentration could clarify this issue. 

In summary, the in-plane thermal conductivity of 
Sr2Cu02Cl2 shows an unusual temperature dependence 
with a pronounced shoulder at high temperatures, similar 
to the behavior found for La2CuC>4 8 and YBa2Cu306- 9 
There is no structural instability in S^CuC^C^. More- 
over, scattering on magnetic excitations is not restricted 
to a narrow temperature interval around the minimum 
of k in these compounds. It is therefore unlikely that 
the double peak structure arises from strong damping of 
the phononic heat current by resonant scattering on soft 
lattice modes or magnetic excitations. The data rather 
indicate a large magnetic contribution to the heat current 
as an intrinsic property of the 2D cuprates. 
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APPENDIX 



To deduce c(T) from the high-temperature series2i 
this series is converted in a series for the entropy s 
in the energy per site e (Ref. l22h . The extrapolations 
are stabilized by information on the ground state en- 
ergy eo = —0.669437(5) (Ref. |23|) . the maximum en- 
tropy s = In 2, and the expected low-energy power law 
s(e) oc (e — eo) 2 ^ 3 . For the latter, Pade approximants are 
applied to s'(e)/s(e) — 2/(3(e— e )) (Dlog Pade approx- 
imation). Very good results are obtained (reliable error 
estimate 10~ 2 from comparing diagonal to non-diagonal 
Pade approximants). The result shown in the inset of 
Fig- El (error 10~ 3 ) is obtained by approximating 

-(e-e y( e )_2l lii / i -eo\ 
s(e) 3J \1 - e J 

which allows for multiplicative logarithmic corrections, 
yielding c(T) cx T 2 (A + hT 7 (l/T)) with A = and 7 = 
1.05(5). We do not, however, exclude a small finite value 
of A as found in spin-wave theory. 
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